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We analyze the difference between the correlation energy as de- 
fined within the conventional quantum chemistry framework and 
its namesake in density-functional theory. Both quantities are rig- 
orously defined concepts; one finds that E^'^' > E^^"^. We give 
numerical and analytical arguments suggesting that the numerical 
difference between the two rigorous quantities is small. Finally, 
approximate density functional correlation energies resulting from 
some popular correlation energy functionals are compared with the 
conventional quantum chemistry values. 



Inroduction 

In quantum chemistry (QC), the exact correlation energy is traditionally de- 
fined as the difference between the exact total energy and the total selfconsistent 
Hartree-Fock (HF) energy: 

rpQC 771 _ rpHF 

-^c,exact ■ ^tot,exact ^tot • \^ ) 

Within the framework of density-functional theory (DFT) ^, on the other 
hand, the correlation energy is a functional of the density E^^'^ [p] . The exact 
DFT correlation energy is then obtained by inserting the exact ground-state 
density of the system considered into the functional E^^'^ [p], i. e. 

rpDFT _ rpDFT r „ 1 /o^ 
^c,exact — [Pexact] ■ \^ ) 



In practice, of course, neither the quantum chemical correlation energy (|I]) nor 
the DFT correlation energy (H) are known exactly. Nevertheless, both quantities 
are rigorously defined concepts. 

The aim of the following section is to give a coherent overwiev of how the cor- 
relation energy is defined in the DFT literature [3-13] and how this quantity is 
related to the conventional QC correlation energy. It will turn out that E^exact 

and -E^f exact generally not identical and that -E^iract ^ ^^e^ct- Furthermore 
we will give an analytical argument indicating that the difference between the 
two exact quantities is very small. 

In the last section we compare the numerical values of approximate conventional 
QC correlation energies with approximate DFT correlation energies resulting from 
some popular DFT correlation energy functionals. It turns out that the difference 
between DFT correlation energies and QC correlation energies is smallest for the 
correlation energy functional of Colle and Salvetti [|14| , [15[| further indicating [|16 



that the results obtained with this functional are closest to the exact ones. 

Basic Formalism 

We are concerned with Coulomb systems described by the Hamiltonian 

H = f + Wcw + V (3) 
where (atomic units are used throughout) 

" ' 1, 



1 ^ 1 

N 

V = ^i;(r,) . (6) 

i=l 

To keep the following derivation as simple as possible, we choose to work with the 
traditional Hohenberg-Kohn fl^ formulation rather than the const rained- search 



representation [T^, |19[ of DFT. In particular, all ground-state wavefunctions 
(interacting as well as non-interacting) are assumed to be non-degenerate. By 
virtue of the Hohenberg-Kohn theorem [T^ the ground-state density p uniquely 
determines the external potential v = v[p] and the ground-state wave function 
\E' [p]. If Vq (r) is a given external potential characterizing a particular physical 
system, the Hohenberg-Kohn total-energy functional is defined as 

E,,[p] = mp]\f+Wcw + Vo\np]) . (7) 

As an immediate consequence of the Rayleigh-Ritz principle, the total-energy 
functional (|^ is minimized by the exact ground-state density Pexact corresponding 
to the potential vq, the minimum value being the exact ground-state energy, i. e. 

Etot,exact -^vq {Pexact\ ■ (8) 



In the context of the Kohn-Sham (KS) scheme the total-energy functional is 
usually written as 



EM=Ts[p]+ J p{r)vo{r) d\ + y / ^ V d V + E., [p] (9) 

where T, [p] is the kinetic-energy functional of non-interacting particles. By virtue 
of the Hohenberg-Kohn theorem, applied to non-interacting systems, the density 
p uniquely determines the single-particle potential Vg [p] and the ground-state 
Slater-determinant 

*''"[P] = ^det {^f/[p]} (10) 
and hence Tg [p] is given by 

Ts[p] = ($^^[p]|T|$^^[p]) 

= E E^f/H(rr(-M^f/[p](r)ciV . (11) 

We mention in passing that the Hohenberg-Kohn theorem can also be formu- 
lated for a "Hartree-Fock world" [^, implying that the HF density uniquely 
determines the external potential. Consequently the HF ground-state determi- 
nant is a functional of the density as well: 

*''^[P] = ^det {^f/[p]} . (12) 

The resulting kinetic-energy functional 

T^^lp] = (<|.^^[p]|f|<|.^^[p]) 




is different from Ts[p\ because the orbitals in (0) come from a /oca/ single-particle 
potential Vs[p] while the orbitals in (|13]) come from the nonlocal HF potential 
v^^[p]. However, the numerical difference between T^^[p] and Ts[p] has been 
found to be rather small [0. 

The remaining term, E^c [p] , on the right hand side of equation (^ is termed the 
exchange- correlation (xc) energy. Comparison of equation (|) with equation (0) 
shows that the xc-energy functional is formally given by 

E.,[p] = {np]\f + Wcitmp])-T4p]-y j^-y^dhd'r' . (14) 

In density-functional theory the exact exchange-energy functional is defined by 
i?r^[p] := ($^^[p]|#c..|$^^[p]) - \ I jP-^0^dhdh' . (15) 



This is identical with the ordinary Fock functional 



b..] = E E 

evaluated, however, with the KS Orbitals, i. e 

E^'^ [p] = E^^ 

The DFT correlation-energy functional is then given by 

E^^^ = E^M-E^'^[p] . 



V^f/ [P] 



(17) 



;i8) 



Inserting the respective definitions ([I^) and (|l^) of E^dp] and E^^'^[p] we find 



^r^p] = mp\\T+wci,\m)-Ts[p\-^ 



P(r)p(r') 



^ff [P] 



(19) 



In terms of the Hartree-Fock total-energy functional 



E^/ b,.] = E E / ^A^r (-^V^) ^,.(r) + / p(r) .;o(r) ciV 



+- 



4,1 j 
1 



ppP^ dh d'r' + ^ [^,.] 



(20) 



and the total-energy functional (0) the DFT correlation energy ([TOD is readily 
expressed as 

Er^[p] = E.M-E;r[^ff{p]] ■ (21) 

By equation (Q), the exact DFT correlation energy is then obtained by inserting 
the exact ground-state density Pexact (corresponding to the external potential Vq) 
into the functional (pT]). By virtue of equation (|[) one obtains 



E 



DFT 

c,exact 



E, 



tot, exact 



E. 



HP 



^1'\P 



exact] 



(22) 



The conventional quantum chemical correlation energy, on the other hand, is 
given by 



E. 



QC 
c, exact 



Etot ,exact 

[Phf] (23) 

where 'pf/'iPHF] are the usual selfconsistent HF orbitals corresponding to the ex- 
ternal potential vo, i. e. phf is that very HF density which uniquely corresponds 
to the external potential vq. Of course, phf and pexact are generally not identical. 
Comparison of ([2^ ) with ( p3D shows that 



E. 



DFT 

c,exact 



E. 



QC 

c, exact 



+ <^ <Mp 



JHF 



E. 



HF 

vo 



^fflp. 



exact 



(24) 



This is the central equation relating the DFT correlation energy to the QC cor- 
relation energy. Since the HF orbitals <^fJ^[pHF] are the ones that minimize the 
HF total-energy functional (EO), the inequality 



E. 



HF 



Vff [phf] 



< E. 



HF 



Vfflp 



exact] 



(25) 



must be satisfied and it follows from equation (|23) that 



^c,exact — ^cexact 



(26) 



Equation ( ^41) tells us that, as a matter of principle, selfconsistent DFT results 
for the correlation energy should not be compared directly with the conven- 
tional quantum chemical correlation energy but rather with the right-hand side of 
equation (p^. In practise, of course, quantum-chemical correlation energies and 
ground-state densities are known only approximately, e. g. , from configuration- 
interaction (CI) calculations. Hence, 



tot,CI 



HF 

tot 



CI\ 



(27) 



is the quantity the selfconsistent DFT correlation energy should in principle be 
compared with. The second term of (27) is readily computed by employing one 
of the standard techniques |]12|, [23, E3|, E3] of calculating the KS potential and its 



orbitals from a given CI density. In the following we shall argue, however, that 
the difference between ^'fexlct £^^^act can be expected to be very small. To 
see this we rewrite equation (|2l) as 



c.exact ^cfixact 



W [Phf] 



E. 



HF 

vo 



3 

'^fa[Px-only] 



E. 



HF 



^ff[Px-anly] 



E. 



HF 

vo 



exact] 



(28) 



where Px-oniy is the ground-state density of an exact exchange-only DFT cal- 
culation ^ and ipf^[px-oniy] are the corresponding KS orbitals. The first 



difference on the right-hand side of equation (|28|) is known to be small p5| , p6[] . 
The second difference, on the other hand, is easily seen to be of second order in 
{px-oniy — Pexact) and is therefore expected to be very small as well: 



E, 



HF 

vo 



^ja [Px-only] 
6E^^ 



HF 



E. 



exact] 



5p(r) 



i^px—onlyij^^ Pexactij^^^ H" 0(ypx—only Pexact) 



Px — only 



Jd^r H- (px-onlyir) - Pexactir)) + 0{px-only " PexactY 
+ 0{px 

—only Pexact) 



The second equality follows from the fact that Px-oniy minimizes the density func- 
tional E^^^ '^ff[p] ■ Hence we conclude that -E^exlrf ~ E^^^^act should be small. 
This estimate is confirmed by results of accurate variational calculations on H~, 
He, Be"*"^, Ne"*"* |T^, as can be seen from Table 1. There, the conventional quan- 
tum chemical correlation energies of these systems are compared with the "exact" 
DFT correlation energies calculated from equation (^2]). For all elements and ions 
shown, the difference between the two quantities is small, as expected. The values 
for H~, He and Be"*"^ also confirm the relation (0). For Ne"*"*^, where the differ- 
ence is only 1 /iH, the inequality (^Bf ) is not satisfied, indicating inaccuracies of 
the numerical procedure involved. 



Table 1: Comparison of DFT correlation energies with conventional quantum 
chemical correlation energies (QC). A denotes the difference between the QC and 
the DFT correlation energy. All numbers in Hartree units. 





DFT 


QC 


A 


H 


-0.041995 


-0.039 821 


+0.002174 


He 


-0.042107 


-0.042 044 


+0.000 063 


Be+2 


-0.044 274 


-0.044 267 


+0.000 007 


Ne+8 


-0.045 692 


-0.045 693 


-0.000 001 



To conclude this section, we mention that there exists yet another possibility of 
defining a density functional for the correlation energy [4-11,13]: 



EM = EM -Elf/ hffip] 



-.HF 



(29) 



where (f^^ [p] are the HF orbitals corresponding to the density p (see equation 
(|I2D). If the exact density Pexact is inserted in (|29| ) ipfflpexact] are the HF orbitals 
corresponding to some unknown external potential Vq whose HF density is Pexact- 
The decomposition 

5o(r) =:i;o(r)+5,(r) (30) 

makes clear that on the single-particle level the definition ( p9D leads to a hybrid 
scheme featuring the ordinary non-local HF exchange potential combined with 
the local correlation potential Vc{tc). In the present paper, this hybrid scheme will 
not be further investigated. 



Correlation Energies from Various DFT Approx- 
imations 

For further analysis, we compare in Tables 0, ^ and |^ the DFT correlation energies 
resulting from various approximations to LYP denotes the correlation- 

energy functional by Lee, Yang and Parr [^, PW91 the generalized gradient ap- 



proximation by Perdew and Wang |^8[, and LDA the conventional local density 
approximation with the parametrisation of E^. by Vosko, Wilk and Nusair [^. 
The first column, denoted by OEP, shows the results of a recently developed 
scheme which employs an optimized effective potential (OEP) including correla- 
tion effects |]TB[. In this scheme the full integral equation of the optimized effective 
potential method [BO, pl[]. 



E/rfV(K'^nr')-w(r'; 

i=i 



I 



E 



J 



^i,{v)ipU^') + c.c. = 
(31) 



with 



:= ;,/f-^rT^ (32) 



is solved semi-analytically by an approved method due to Krieger, Li and lafrate [32 



33, 34 



1 . . - 
K'^J'^ir) ^ Kt^'(r) = -7-T E P^^(r) h«.(r) + {K^'ff - n..,.)] (33) 

P'^K'^) i=l 

where the constants {v^J"/' — Uxda) are the solutions of the set of linear equations 
E (5,. - M,,.) iy^l'' - n.,,.) = y^., - u,,,^ J = 1, . . . , AT, (34) 



=1 

with 



Pja(r)pia(r 



(36) 



i=l P(t{'^) 



Here, Uxcja denotes the average value of Uxcjo-(i') taken over the density of the ja 
orbital, i. e. 

Uxcja = J pjaiT^)uxcja{r)d^r (37) 

and similarly for V^cja- Like in the conventional Kohn-Sham method, the xc- 
potential resulting from equation (p3D leads to a single-particle Schrodinger equa- 
tion with a local effective potential 

(-^ + ^o(r) + / rfV + V;^r (r) j ^,.(r) = 5,.y.,.(r) (38) 

The selfconsistent solutions v^jo-(r) of equation ( pH]) with lowest single-particle 
energies ej^^ minimize the total-energy functional 

E^^'i^,.] = E E /^L(r) (-M ^..(r)rf^- 



+ y p(r) 'yo(r) (i^r 



1 /■ /■ p(r)p(r') 3 3 , 



|r — r 



-i E E lld'rd 

■^CS 



^3^ ^3^/ V^jV('^)y^L(r0^fca(r)¥^ja(rO 



+i^c"' K^.J] • (39) 



In the above equation, E^^ denotes the CoUe-Salvetti functional JT^ |T5| for the 
correlation-energy given by 



= -a6|7(r)e(r) 



where 



7(r) 
r/(r) 



4 PrWPiW 

p(r)2 ^ 
1 + dp{r)~3 

P(r) 



'36" 



'Cp(r) 3 



7]{r) 



(40) 

(41) 
(42) 
(43) 



The constants a, b, c and d are given by 

a = 0.04918, 6 = 0.132, 
c = 0.2533, d = 0.349. 

In Table ^, the four approximate DFT correlation energy functionals are evaluated 
at the exact densities [O] of H~, He, Be+^, Ne+^ and compared with the exact 



DFT correlation energies given by equation (^2]). On average, the OEP values 
are superior. 

In Table ^ selfconsistent DFT correlation energies are compared with QC values 
taken from . In these selfconsistent calculations the approximate correlation- 
energy functionals E^^^, E^^^^, E^^"^ are complemented with the approximate 
exchange-energy functionals E^^^ [^, E^^^^ |2^ and E^^^, respectively. In the 
OEP case, the DFT exchange-energy functional {^i7\) is of course treated exactly. 
The numerical data show three main features. 

1. For most atoms, the absolute value of Ef^ is smaller than the absolute cor- 
relation energy obtained with any DFT method, as it should be according 
to the relation (^). 

2. The values of E^^^, E^^^ , j^PWdi j^qc g^gj-gg quite closely with each 
other while the absolute value of E^^^ is too large roughly by a factor of 
two. We mention that due to the well known error cancellation between 
E^^^ and E^^^, the resulting LDA values for total xc energies are much 
better. 

3. The difference between E^^'^ and E^^ is smallest for the E^^^ values, 
larger for E^^^ and largest for E^^^^. The difference between E^*^' and 
fjDFT j^g^g three sources: 



Table 2: Non-relativistic absolute correlation energies resulting from various 
approximate DFT correlation energy functionals, evaluated at the exact ground- 
state densities of the respective atoms. \A\ denotes the mean absolute deviation 
from the exact DFT correlation energy. All numbers in Hartree units. 





OEP 


LYP 


PW91 


LDA 


EXACT 


H" 




0.0297 


0.0299 


0.0320 


0.0718 


0.0420 


He 




0.0416 


0.0438 


0.0457 


0.1128 


0.0421 




0.0442 


0.0491 


0.0537 


0.1512 


0.0443 


Ne 


+8 


0.0406 


0.0502 


0.0617 


0.2030 


0.0457 




A 




0.0045 


0.0058 


0.0097 


0.0912 





Table 3: Non-relativistic absolute correlation energies of first and second row 

atoms from selfconsistent calculations with various DFT approximations. QC 
denotes the conventional quantum chemistry value. | A|% denotes the mean value 



of 



DFT 

c 



E^'") / E'^^ in percent. All other numbers in Hartree units. 





OEP 


LYP 


PW91 


LDA 


QC 


He 


0.0416 


0.0437 


0.0450 


0.1115 


0.0420 


Li 


0.0509 


0.0541 


0.0571 


0.1508 


0.0453 


Be 


0.0934 


0.0954 


0.0942 


0.2244 


0.0943 


B 


0.1289 


0.1287 


0.1270 


0.2906 


0.1249 


C 


0.1608 


0.1614 


0.1614 


0.3587 


0.1564 


N 


0.1879 


0.1925 


0.1968 


0.4280 


0.1883 





0.2605 


0.2640 


0.2587 


0.5363 


0.2579 




P 


0.3218 


0.3256 


0.3193 


0.6409 


0.3245 


Ne 


0.3757 


0.3831 


0.3784 


0.7434 


0.3905 


Na 


0.4005 


0.4097 


0.4040 


0.8041 


0.3956 


Mg 


0.4523 


0.4611 


0.4486 


0.8914 


0.4383 


Al 


0.4905 


0.4979 


0.4891 


0.9661 


0.4696 


Si 


0.5265 


0.5334 


0.5322 


1.0418 


0.5050 




P 


0.5594 


0.5676 


0.5762 


1.1181 


0.5403 




3 


0.6287 


0.6358 


0.6413 


1.2259 


0.6048 


CI 


0.6890 


0.6955 


0.7055 


1.3289 


0.6660 


Ar 


0.7435 


0.7515 


0.7687 


1.4296 


0.7223 




A 


% 


3.13 


4.52 


5.10 


120 





Table 4: Non-relativistic absolute correlation energies of atoms from selfcon- 
sistent calculations with various DFT approximations. All numbers in Hartree 
units. 





OEP 


LYP 


PW91 




OEP 


LYP 


PW91 


K 





8030 





7821 





7994 


Rb 


1.7688 


1.7832 


1.9509 


Ca 





8269 





8329 





8467 


Sr 


1.8222 


1.8355 


2.0056 


Sc 





8832 





8855 





9033 


Y 


1.8763 


1.8863 


2.0671 


Ti 





9371 





9374 





9613 


Zr 


1.9281 


1.9363 


2.1307 


V 





9882 





9882 


1 


0198 


Nb 


1.9475 


1.9558 


2.1899 


Cr 


1 


0073 


1 


0086 


1 


0736 


Mo 


1.9905 


2.0003 


2.2551 


Mn 


1 


0812 


1 


0861 


1 


1375 


Tc 


2.0796 


2.0874 


2.3412 


Fe 


1 


1597 


1 


1620 


1 


2158 


Ru 


2.1571 


2.1637 


2.4254 


Co 


1 


2324 


1 


2331 


1 


2933 


Rh 


2.2278 


2.2340 


2.5081 


Ni 


1 


3009 


1 


3010 


1 


3700 


Pd 


2.3123 


2.3154 


2.6074 


Cu 


1 


3693 


1 


3694 


1 


4562 


Ag 


2.3561 


2.3649 


2.6705 


Zn 


1 


4273 


1 


4303 


1 


5212 


Cd 


2.4146 


2.4247 


2.7373 


Gel 


1 


4704 


1 


4753 


1 


5768 


In 


2.4600 


2.4704 


2.7964 


Ge 


1 


5101 


1 


5174 


1 


6343 


Sn 


2.5024 


2.5135 


2.8577 


As 


1 


5465 


1 


5570 


1 


6917 


Sb 


2.5419 


2.5544 


2.9193 


Se 


1 


6177 


1 


6288 


1 


7662 


Te 


2.6134 


2.6252 


2.9965 


Br 


1 


6795 


1 


6912 


1 


8393 


I 


2.6763 


2.6876 


3.0726 


Kr 


1 


7355 


1 


7493 


1 


9112 


Xe 


2.7338 


2.7456 


3.1475 



a) The values of E^'^ are only approximate, i. e. not identical with 



■'-'c, exact- 



(b) The values of are only approximate, i. e. not identical with 

pDFT 
c, exact' 

(c) As shown in the last section, the exact values E^^^^^ct ^iid -Efexoct 
not identical. 

Currently it is not known with certainty which effect gives the largest con- 
tribution. However, with the arguments given in the last section, we expect 
the contribution of (c) to be small. Assuming that the quoted values of 
E^'-^ are very close to Ef exact we conclude that E^^^ is closest to -Ef exact- 

Table ^ shows correlation energies of atoms K through Xe obtained with the 
various selfconsistent DFT approaches. In almost all cases, the absolute OEP 
values for E^ are smallest and the ones from PW91 are largest, while the LYP 
values lie in between. In most cases, E^^^ and E^"^^ agree within less than 1 % 
while is larger (by up to 10 %) as the atomic number Z increases. We 

emphasize that reliable values for E^'-^ do not exist for these atoms. 
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